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Abstract 

We study the energy levels of H2 molecules in a superstrong magnetic 
field {B ^ 10^^ G), typically found on the surfaces of neutron stars. The 
interatomic interaction potentials are calculated by a Hartree-Fock method 
with multi-configurations assuming electrons are in the ground Landau state. 
Both the aligned configurations and arbitrary orientations of the molecular 
axis with respect to the magnetic field axis are considered. Different types of 
molecular excitations are then studied: electronic excitations, aligned (along 
the magnetic axis) vibrational excitations, transverse vibrational excitations 
(a constrained rotation of the molecular axis around the magnetic field line). 
Similar results for the molecular ion are also obtained and compared with 
previous variational calculations. Both numerical results and analytical fitting 
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formulae are given for a wide range of field strengths. In contrast to the zero- 
field case, it is found that the transverse vibrational excitation energies can be 
larger than the aligned vibration excitation, and they both can be comparable 
or larger than the electronic excitations. For B ^ Bcrit = 4.23 x 10^^ G, 
the Landau energy of proton is appreciable and there is some controversy 
regarding the dissociation energy of H2. We show that H2 is bound even for 
B » Bcrit and that neither proton has a Landau excitation in the ground 
molecular state. 

32.60.+i, 97.10.Ld, 3L15.+Z, 97.60.Jd 
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I. INTRODUCTION 



Since the pioneering work of Schiff and Snyder [|T|, especially during the last 20 years, 
there has been considerable interest in the properties of matter in a strong magnetic field. 
While the early studies were mainly motivated by the fact that high magnetic field con- 
ditions can be mimicked in some semiconductors where a small effective mass and a large 
dielectric constant reduce the electric force relative to the magnetic force, the recent interest 
in this problem has been motivated by the huge magnetic field ~ 10^^ G already discovered 
in many neutron stars and the tentative suggestion for fields as strong as 10^^ G. The sur- 
face layer of these neutron stars then consists of highly-magnetized matter. Understanding 
the physical properties of atoms, molecular chains, and condensed matter in fields of such 
extreme magnitude (see Ref. [§] for an early general review and @] for a recent text on atoms 
in strong magnetic fields) is important for interpreting the radiation from the neutron stars 
that may be observed in the present and future X-ray satellites (e.g., 0), therefore provides 
important information about the internal structure of neutron stars. 

In superstrong magnetic fields the structure of atoms and condensed matter is dramati- 
cally changed by the fact that the magnetic force on an electron is stronger than the Coulomb 
force it experiences, i.e., the electron cyclotron energy (the Landau energy level spacing) 

eB 

hoje = h = 11.57Bi2 keV, (1.1) 

THeC 

where B12 is the magnetic field strength in units of 10^^ G, is much larger than the typical 
Coulomb energy. In the direction perpendicular to the field, the electrons are confined to 
move on cylindrical Landau orbitals around a nucleus. The orbitals have radii 

Pm = i2m+iy/^p, m = 0,l,2,---, (1.2) 

where p is the cyclotron radius 

P = = ao (f )'^' = 2.57 X 10-'' B;,'/' cm. (1.3) 
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Here = Ti^ / {m(,c^) is the Bohr radius and Bo is the atomic unit for the magnetic field 
strength, 

Bo = = 2.35 X 10^ G, b = ^ = 425 5i2. (1.4) 

h Bo 

Throughout this paper we consider strong fields in the sense of 6 >> 1, so that the Coulomb 
forces act as a perturbation to the magnetic forces on the electrons, and the electrons are 
confined to the ground Landau level (so called "adiabatic approximation" [0]). Because of 
this extreme confinement of electrons in the transverse direction, the Coulomb force becomes 
much more effective for binding electrons in the parallel direction, therefore giving greatly 
increased binding energy. The atom has a cigar-like structure. Moreover, it is possible for 
these elongated atoms to form molecular chains by covalent bonding along the field direction 

ii- 

Significant efforts have been devoted to the theoretical study of atoms in a superstrong 
magnetic field (^ 10^^ G) [Q. The methods that have been employed include variational 
calculations(e.g., |^), Thomas- Fermi-type statistical models P, density functional theory 



0, and self-consistent Hartree-Fock method p0|-p!2|, which we consider to be the more 
theoretically justified and reliable a method. Accurate calculations of the energy levels of 
the H atom in magnetic fields of arbitrary strength have also been performed [ll^. By 
contrast, there are only limited studies on molecules in superstrong magnetic field; nearly 
all of these focus on the molecular ion ( [p!^ -[T7| and references therein). As is 
unstable against forming H2, understanding the physical properties of H2 molecule is of 
greater practical interest, since H2 is likely to exist in the atmosphere of sufficiently cool 
neutron stars ||T8| , p!9| 

We have recently calculated the ground state binding energies of different forms of hy- 
drogen (H, H~, H^, H2, H3, ■ ■ -, Hoo) in a strong magnetic field B ^ 10^^ G ( [|], hereafter 
referred as Paper I). In particular, for the first time, reliable electronic dissociation energy 
of H2 molecule in magnetic field of such magnitude was obtained. In this paper, we extend 
our study to consider various excitation levels of the molecule. 



In the zero-field case, to study the molecular spectra, one usually uses Born-Oppenheimer 
approximation to separate the motion of the ions from that of the electrons. Such a procedure 
is valid if the electronic energy-level spacings are large compared to the typical energy-level 
spacings associated with the ion motion. In a strong magnetic field, however, the separation 
of motion becomes much more complicated, even for the hydrogen atom p]|-p2[]. Moreover, 



as we shall see, in a superstrong magnetic field, the energy-level spacings associated with the 
vibrations of the ions can be comparable to or even larger than the spacings of the electronic 
excitations. In this paper, we will use the standard Born-Oppenheimer approximation and 
focus on calculating the interatomic interaction potential for fixed ion positions (Sec. III). 
We then obtain the molecular excitation levels based on this potential curve (Sec. IV). As in 



the case of a neutral atom pO| , it is convenient to define a critical field strength by equating 
the cyclotron energy of the proton hup = h^eB/nipc) to the typical electronic excitation 
energy (~ In 6 in atomic units), i.e., 

bcrit = — \nbcrit = 1.80 X 10^ B,rit = KritBo = 4.23 X 10^=^ G. (1.5) 

rrie 

We shall give quantitative results for the regime Bo « B « Bcrit in Sec. II- IV, using 
the standard Born-Oppenheimer procedure. Rigorous calculations for the molecule when 
B ^ Bcrit, taking account of the quantum mechanics of the proton motion, are difficult. 
Nevertheless, in Sec. V we shall describe an approximate solution to the four-body problem 
of the H2 molecule in the B » B^rit regime, where the effects of finite proton mass on the 
electronic states and the energies of the molecule are strong, and we give a rigorous lower 
limit to the ground-state dissociation energy. 

Throughout this paper, we shall use nonrelativistic quantum mechanics, even for ex- 
tremely strong magnetic field, B > Brei = [hc/e^YBo = 4.414 x 10^^ G (note that Brei is 
close to Bcrit only by coincidence), at which the transverse motion of the electron becomes 
relativistic. The nonrelativistic treatment of bound states is valid for two reasons: (i) The 
energy of a relativistic free electron 



E 



cY, + my (1 + 2nLB/Brei) , (1.6) 



where pz is the hnear momentum along the field axis, n^, is the quantum number for the 
Landau excitations, reduces to ^ meC^ +p^/(2me) as along as the electron remains in 
the ground Landau level and nonrelativistic in the 2;-direction; (ii) The shape of the Landau 
wavefunction in the relativistic theory is the same as in the nonrelativistic theory (as we 
see from the fact that p is independent of mass). Therefore, as long as Eb/ {Tnf,(?) « 1, 
where Eb is the binding energy of the bound state, the relativistic effect remains a small 



correction [p3 



The paper is organized as follows. In Sec. II we consider some general features and 
approximate scaling relations for various excitation levels. Sec. Ill contains a detailed de- 
scription of our method for calculating the interatomic interaction potential. The numerical 
results and fitting formulae for the molecular excitation levels are presented in Sec. IV. We 
study the electronic structure of the molecule in the B » Bcru regime and consider the 
effects of finite proton mass on the energies in Sec. V. Our general conclusion is presented 
in Sec. VI. Appendix A summarizes some useful mathematical relations for the Coulomb 
integrals of Landau functions, and in Appendix B we discuss a refined method for calculating 
the electronic energy of ^2 for general orientation of the molecular axis. 



II. QUALITATIVE DISCUSSION AND APPROXIMATE SCALING RELATIONS 

FOR EXCITATION ENERGIES 

In a superstrong magnetic field satisfying 6 >> 1, the spectra of a single H atom can be 
specified by two quantum numbers (m, u), where m measures the mean transverse distance 
(Eq. [1.2]) of the electron to the proton, while u is the number of nodes of electron's z- 
wavefunction (along the field direction). The wavefunction of the (m, u) state in cylindrical 
coordinates (p, 0, z) is given by 

^mu = W^{r^)U,{z), (2.1) 
where Wm is the ground state Landau wavefunction 
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W^{v^) ^ Wm{pA) = -^L=(^Ye-^'/^P"e-'^t (2.2) 



The states with z/ 7^ resemble a zero-field hydrogen atom with small binding energy 
\Ey\ ^ l/(2z/2) and we shall mostly focus on the tightly-bound states with z/ = 0. For 



the ground state (0,0), the sizes and of the atomic wavefunction perpendicular and 
parallel to the field and the binding energy \Ea\ (in atomic units) are given by 

^ ^ ^ 0.16/2; / = ln6. (2.3) 



For the tightly-bound excited states (m, 0) we have similar relations but with p replaced by 
[(2m + l)/&]^/2 and / replaced by /,„ = ln[6/(2m + 1)], so that 

E,(m)~-0.16C (2.4) 

Recall that in atomic units, mass is in units of the electron mass me, energy is expressed in 
units of e^/ao = 2 Ry, length is in units of Bohr radius oq, and the units for magnetic field 
is Bo (Eq. [1.4]). The numerical factor 0.16 in Eqs. (2.3)-(2.4) is an approximate value for 
B12 ^ 1. For convenience, accurate numerical results for Eaim) are listed in Table I. [| 

In a superstrong magnetic field, the mechanism of forming molecules is quite different 
from the zero- field case (Paper I, [Q). The spins of the electrons of the atoms in a strong 
magnetic field are all aligned anti-parallel to the magnetic field, and therefore two atoms in 
their ground states do not easily bind together according to the exclusion principle. Thus 
two H atoms, both in the m = ground state, do not form tightly bound molecule. Instead, 
one H atom has to be excited to the m = 1 state. The two H atoms, one in the ground state 
(m = 0), another in the m = 1 state then form the ground state of H2 molecule by covalent 



more accurate fitting formula for the ground state binding energy of H atom is \Ea\ = 0.16 A l'^, 

with 

1 + 1.36 X 10-2 [ln(1000/6)]2-5, if 6 < 10^; 

A ; 

1 + 1.07 X 10-2 [ln(Vl000)]i-6, if 6 > 10^ 
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bonding. Since the "activation energy" for exciting an electron in the H atom from Landau 
orbital m to (m + 1) is small (see Eq. [2.4]), the resulting molecule is stable. The interatomic 
separation Zg and the dissociation energy D of the H2 molecule scale approximately as 

^o = e^.~y, (2-5) 

where the dimensionless factor ^ decreases very slowly with increasing B (e.g., C, — 2.0 for 
B12 = 0.1 and ^ ~ 0.75 for Bu = 100; see Table I of Paper I and our Eq. [5.2]). 

Another mechanism of forming H2 molecule in a superstrong magnetic field is to let both 
electrons occupy the same m = Landau state, while one of them occupies the z/ = 
orbital and another u = 1 orbital. This costs no "activation energy". However, the resulting 
molecule tends to have small dissociation energy, of order a Rydberg. We shall refer to 
this electronic state of the molecule as the weakly-bound state, to the states formed by two 
electrons in the u = orbitals as the tightly-bound states. As we will see below, as long as 
/ >> 1, the weakly-bound state only constitutes an excited energy level of the molecule. 

We now consider various molecular excitations and derive approximate scaling relations 
for the excitation energies. 



A. Electronic Excitations 

The electronic excitations of are similar to those of the H atom, namely the electron 
can occupy different m Landau orbitals. Thus m = is the ground state, m = 1, 2, ■ ■ ■ are 
the excited states (although they are not necessarily bound relative to the free atom in the 
ground state). 



^In several recent papers [^] on the molecular binding in strong magnetic field, Korolev and 
Liberman failed to identify the tightly-bound states. Also, their variational calculation of the 
weakly-bound state significantly underestimates the binding energy because it neglects the over- 
lapping of the electron wavefunctions. As a result, their claim that hydrogenlike gas in strong 



magnetic field can form Bose-Einstein condensate is incorrect (see also [26,^7]]). 



There are two types of electronic excitations in H2. (i) The electrons can occupy different 
orbitals other than the ground state {mi, 1712) = (0,1), giving rise to the tightly-bound 
{u = 0) electronic excitations. For example, the first excited level is (0,2), the second 
excited level is (0,3), etc. The number of single m-excitation states (mi, 7712) = (0, m2) 
which are bound relative to two isolated H atoms in the ground state is expected to increase 
as the magnetic field increases. Double m-excitations are also possible, but as we shall see, 
they are bound only when the magnetic field strength is much higher than 10^^ G. The 
energy spacing between the two adjacent electronic states (0,m) and (0,m + 1) is 



Thus as m increases, the energy spacing decreases, (ii) The molecule is formed by two 
electrons in the (m, u) = (0, 0) and (0, 1) orbitals. The dissociation energy of this weakly- 
bound state is of order a Rydberg, and does not depend sensitively on the magnetic field 
strength. Note that for relatively small magnetic field {B12 ^ 0.2), the weakly-bound state 
actually has lower energy than the tightly-bound states (see Sec. IV. A), i.e., b > 10^ is 
required for the "strong field" regime to apply fully. 



In the Born-Oppenheimer approximation, the motion of the two protons is governed by 
the interatomic potential U{Z,R±), i.e., the electronic energy when the relative positions 
of the protons are kept at Z along the field direction and R± perpendicular to it. We 
first consider the aligned vibrational excitations for oscillations of Z about the equilibrium 
separation Zg. For this purpose we need to estimate the excess potential SU{6Z) = U{Zo + 



Since Zg is the equilibrium position, the sum of the first order terms in 6Z, coming from 
proton-proton, electron-electron, proton-electron Coulomb energies and quantum mechanical 
electron kinetic energy, must cancel. Thus we have 6U oc {5Zy for small 5Z. Consider 




(2.6) 



B. Aligned Vibrations 



SZ, 0)-U{Zo,0). 
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various contributions to the energy of the molecule: The proton-proton interaction is l/Z 
(in atomic units) without a logarithmic factor; but the dominant contribution is the proton- 
electron Coulomb energy ~ //Z, where the logarithmic factor / >> 1 comes from the 
Coulomb integral over the "cigar-shaped" electron distribution. Both / and change as 

+ bZ, but the largest change comes from the quadratic term ~ 
Thus the excess potential is of order 

5f/(5Z) ~ / ^ ~ (r^/^) {bZf. (2.7) 

In atomic (electron) units the reduced mass of the proton-pair in H2 is /i = mp/(2me), 
where mp and me are proton and electron mass (for HD the factor 1/2 is replaced by 2/3). 
For small-amplitude oscillations in the potential of equation (2.7), we obtain a harmonic 
oscillation spectrum with excitation energy quanta Tiios^ given by 

~r'/2^V"'/', (2.8) 

for a molecule in the ground electronic state. The scaling with B of Tiuj^]^ is thus almost the 
same as the dissociation energy D in Eq. (2.5). The number of aligned vibrational levels is 
'n\\max ~ D /huj\\ ~ (^/^)^^^) where ^ decreases even more slowly with increasing field strength 
than does. 

C. Transverse Vibrations 

The strong magnetic field breaks the rotational symmetry for the molecular axis and, 
instead of rotations of the field-free case we have oscillations in the two-dimensional plane 
of the R_L vector^. The degeneracy in the azimuthal angle is still retained. To study the 



Strictly speaking, the transverse vibration and the aligned vibration are coupled, and they are 
governed by the two dimensional potential U{Z,R±). Since the transverse vibrational excitation 
energy is larger than the aligned vibrational excitation, the timescale for the protons to adjust their 
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transverse vibration spectrum, we need to estimate the order of magnitude of the excess 
potential 5U{R±) = U{Zo, R±) - U{Zo, 0). 

As mentioned before, the factor / in the expression l/Zg for the dissociation energy D 
(Eq. [2.5]) comes from a Coulomb integral over the electron charge distribution. This integral 
is of the form \r{Lz/p), where p — b"^^'^ is the typical size of the electron wavefunction 
perpendicular to the field for R± = 0. When the protons are displaced by R± from the 
electron distribution axis, the Coulomb integral can be approximately obtained by replacing 
p with (p^ + i?^)^/^. Our order of magnitude expression for 5U is then 

SU{Ri_) ~ In (l + r'^i) ~ C^l In (l + 6 RI) . (2.9) 

Equation (2.9) holds for any R± « Zo ~ but it can be approximated by a quadratic 
expression for the small- amplitude case of R± ^ p — << Zg. In this approximation 

we have SU ~ ^~^lbR']^. The energy quanta for the small- amphtude transverse vibration is 
then 

^^±0- (r'/^^)'^%"'/', (2.10) 

where the subscript indicates that we are at the moment neglecting the magnetic forces 
on the protons which, in the absence of Coulomb forces, lead to the cyclotron motions of 
the protons. Note that huj±o in Eq. (2.10) increases as 6^^^ with increasing field strength, 
faster than the logarithmic behavior of hu!\\ and D, but slower than the linear behavior 
of the cyclotron energy. For sufficiently large 6 >> 1 we have tiuj±o » huj\\. However, 
the quadratic harmonic oscillator approximation is valid only for R'j_ up to ~ = 
i.e., for 5U only up to SUho ~ which is less than the maximum possible potential 

^Umax ~ D ~ The number of harmonic oscillation levels in the quadratic regime is 

then 



positions is much longer than the timescale for oscillations with R± ^ and we can consider 
transverse vibrations with fixed values of Z. However, since 5Z « Z, an approximate separation 
is possible with Z replaced by Zg for the transverse vibrations. 
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n±ho ^ 




1/2 



(2.11) 



The degeneracy of the n^-th harmonic oscillation level is n±. For n±ho >> 1, the statistical 
weight of all harmonic oscillation levels is of order {n±ho)'^- If we neglect the difference 
between ^ and unity (and between n and mp/me), we see that n±ho would be less than unity 
when B ^ Bcru, where Bcru is defined in Eq. (1.5). 

We now consider large amplitude transverse oscillations assuming that the magnetic 
force on the proton can be neglected. For a transverse oscillation wavef unction where the 
maximum value Rmax of R± (the outer classical turning point) satisfies p ^ Rmax ^ Zo, 
we must use the logarithmic form of Eq. (2.9) for the potential SU{R±). The energy level- 
spacing decreases with increasing Rmax- We can calculate the number of nodes n±{Rmax) 
of the wavcfunction as a function of Rmax from a WKB integral of the wave number k{R±) 
over dR±. Since we only need an order of magnitude estimate, we replace the integral by 



The maximum number of nodes n±max can be obtained by substituting Zq ~ ^/Z for Rmax- 
Neglecting In / compared with / itself, we have n±max ~ (^A*)^^^; independent of field strength 
and the same order of magnitude as n\\max- 

Because of the azimuthal symmetry in the two-dimensional R^-plane, the total statistical 
weight of the transverse excitation levels is ~ ""•imaa; ^ ^1^- If ^ >> bcrit, n±ho in Eq. (2.11) 
would be much less than unity and the zero-point energy e±zp, i.e., the spacing between the 
lowest levels, is not given by Eq. (2.10). Formally, one could use Eq. (2.9) and estimate 
the zero-point vibration amplitude as the value of Rmax for which Eq. (2.12) gives n± = 1. 
This would give a zero-point energy which is less than D, but this expression is incorrect, 
since the neglect of the magnetic forces on the protons is unjustified when B » Bcru- The 
cyclotron energy of the proton is hUp = heB / {nipc) = {rrie/mp) b (a.u.). The ratio ujp/uj±Q is 
of order (hrrie/lmpY/'^ — {h/hcritY^"^ (omitting the factor ^). When Uujp is much larger than 



k{Rmax)Rmax, where k{R±) ~ [fxSU{R±)Y/^. Using Eq. (2.9) this gives 



r ~\ 1/ 2 

n^{Rmax)^ fiC'lH^ + bRLa.) R, 



■max- 



(2.12) 
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hu}±o, the magnetic forces on the protons are important. We will return to this subtle issue 
in Sec. V. 



III. METHODS FOR CALCULATING THE INTERATOMIC POTENTIAL 

In the Born-Oppenheimer approximation, the interatomic potential U{Z, R±) is given by 
the total electronic energy E{Z, R±) of the system when the relative positions of the protons 
are Z along the field direction and R± perpendicular to it. Once E{Z, R±) is obtained, the 
electronic equilibrium state can also be determined by locating the minimum of the E{Z, 0) 
curve. 



A. The Aligned Case: R± = 

Our method for calculating E{Z, 0) is the same as in Paper I. It can also be used to 
obtain the energy curves for the excited electronic states. Here we summarize and extend 
our method to take account of "configuration interaction" in H2 more accurately. 



1. H2 Molecular Ion 



For 112", the Hamiltonian for the electron is 

92 

Ho^Hb-- — -^-^ , (3.1) 

zme OZ"^ Ta Tb 

where ta and tb are the distances between the electron and the two fixed protons, located 
Ski z = ±Z/2 along 2;-axis. In Eq. (3.1), Hb is the magnetic part of the Hamiltonian 

= 77^ f P± + - A y + • S, (3.2) 

zme V c / mgC 

where A = B x r/2 and S is the electron spin operator. Note that for electrons in the 
ground Landau level, we have 

HB[Wm{v^)x{i)] = 0, (3.3) 
13 



where xii) is the electron spinor with the spin ahgned in the — z-direction (anti-parallel 
to the field). Thus we can set Hb — 0. With the electron wavefunction given by 
^mo(r) = W^m(r±)/mo(2;), we average over the transverse direction and obtain a one- 
dimensional Schrodinger equation 

fmO -Vrn{z)fmO — ^mofmO- (3-4) 



Here the averaged potential is given by 



Vm{z)^ ydV|W^m(r±)|'(^ + ^)=V;n(^-|)+Kn(^+|), (3.5) 

where 

Vm{z) = /dV|W^m(r±)r^, (3.6) 

which can be evaluated numerically (Paper I). In Eqs. (3.4)-(3.6) and hereforth we employ p 
as the length unit in all wavefunctions and average potentials (except otherwise noted) . We 
solve the eigenvalue Smo by integrating Eq. (3.4) numerically from z = +00 to ^ = subject 
to appropriate boundary conditions (Paper I). The total electronic energy is then given by 

E{Z,0)^emo + ^. (3.7) 

Clearly, m = is the ground state, while m = 1, 2, • • • are the excited electronic states. 

We also note that the excited state of H2" in which the electron occupies the u > orbital 
is not bound relative to the free atom in the ground state. 

2. H2 Molecule: Tightly-Bound States (m, v) = (mi , 0) , (m2 , 0) 

For H2, we use the Hartree-Fock (HF) method to take account of the interaction between 
the electrons. The Hamiltonian of the system is 

H = Ho{l) + Ho{2) + — + ^, (3.8) 

ri2 Z 
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where Ho is given by Eq. (3.1) and ri2 = |ri — r2|. For the (777,1,777.2) electronic state 
{mi 7^ 7772), the two basis wavefunctions (orbitals) for the electrons are 

^mio(r) = WmAr^)fm,Q{z), (3.9) 

$m.o(r) = W^,{v,_)f^,^{z). (3.10) 
The two-electron wavefunction is then given by 

*(ri, rs) = ;^[$mio(ri)$m2o(r2) - ^mA^2)^mA^l)]- (3.11) 

After averaging over the transverse direction, the standard HF equations reduce to a set of 
one-dimensional equations for /^^q ^-nd /m20- 



Vm{z) —Km{z) - Sr. 



o2 
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/mo(^) = -T--'m(^), 777 = 7771,7772, (3.12) 



2mefP' dz'^ p p 

where is given by Eq. (3.5); the direct and exchange potentials K and J are given by 

KmA^) = j dz' f^^Q{z'fD^^rn2{z - z'), (3.13) 
Jmi{z) — fm2o{z) J dz' fmio{z')fm2o{z')Ejnim2 

{z-z'), (3.14) 

and similarly for K^^ and Jm2- In Eqs. (3.13)-(3.14), D^^^^ and £^mim2 the direct and 
exchange interaction kernels defined by 

Dm,m2{zi -Z2)^ [ 0?V±0?W | W^mi (ri±) | (r2±) l' — , (3.15) 
J 7-12 

Em,m2{zi-Z2)= /o? WW^^mi (ri±)Pr„2 (r2±) W^^, (r2±) P^^, (rix) — ■ (3.16) 

7*12 

The functions Dmj^rn2{z) and ii^mim2(-2) are related to the Coulomb interaction potential 
(Eq. [3.6]) by 



mi+m2 



D^^rn2{z)= ^^*^"^^'"^2);^V;(^--^^, (3.17) 



mi+m2 



1 



^mim2(^) = H ^'*^"^^'"^^'*72^'(^)' 

where the coefficients dg and are given in Paper I. We solve Eq. (3.12) numerically us- 
ing a shooting algorithm (for detail, see Paper I). Once the wavefunction fmo{z) and the 
eigenvalues Emo are obtained, the total electronic energy of the system is calculated via 

15 



E = {^\H\<if) 

^ /" 2 2 

= — + £miO + ^m20 T dZidZ2fmio{zi) /m2o(^2) -Dmim2(^l ~ ^2) 

Z p J 

+ — dZidZ2fm io(^l)/m2o(^2)/mio(^2)/m2o(^l)-^mim2 (^1 ^2)1 (3.19) 
P -I 

where the 4th term on the right hand side represents the electron direct interaction (— i?'^*'"), 
and the 5th term the exchange interaction [—E'^^'^^). 

The Hartree-Fock method discussed above can be used to obtain accurately the electronic 
energy near the equilibrium separation Zo- However, as noted in Paper I, as Z increases, the 
resulting E{Z, 0) becomes less reliable. Moreover, as Z ^ oo, E{Z, 0) does not approach 
the sum of the energies of two isolated atoms, one in the mith state, another in the m2th 
state. The reason is that as Z increases, a second configuration of electron orbitals becomes 
more and more degenerate with the first configuration in Eq. (3.11), and there must be 
mixing of these two different configurations. This "configuration interaction" also occurs 
in the zero-field H2 molecule [^. Here the electron configuration that mixes with \E'i = \& 



(Eq. [3.11]) is 

^2(ri, ra) = -^[<l>™,i(ri)$^,i(r2) - $™,i(r2)$m2i(ri)], (3.20) 

which is the same as \E'i except v = 1 m the electron orbitals. Both \E'i and ^2 have the 
same symmetry with respect to the Hamiltonian in Eq. [3.8]: the total angular momentum 
along the 2;-axis is M^^ = 1, the total electron spin is Msz = —1, and both \E'i and \E'2 are 
even with respect to the operation —>■ — r^. As a result, the matrix element |\1'2) is 

nonzero. 

To take account of the mixing of these two configurations, we need to extend the standard 
HP method involving one configuration to HF with multi- configurations (HPMC). This is 
done as follows. We calculate the energies and wavefunctions for both and \E'2 using the 
HP equations (Eq. [3.12] with v = Q and z/ = 1). The matrix elements Hij {i,j = 1,2) 
are then calculated: Hu = (\E'i|if|^i) is given by Eq. (3.19) and the expression for H22 = 
(\E'2|-f^|^2) is similar. The mixing matrix element is given by 
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//l2 = (*l|//|*2) 




— /(i^l(i^2/mio(^l)/mil(^l)/m2o(^2)/m2l(^2)-Dmim2(^l ~ ^2) 




^17111712 

[Zl - Z2 . 

P J 



(3.21) 



The total electronic energy in this HFMC scheme is obtained by solving the secular equation 
det \Hij — E5ij\ = 0, which yields, for the lowest energy state 



In Fig. 1, we show the electronic energy curves of H2 at Bu = 1, obtained using our 
HFMC method. The tightly-bound electronic states are (mi, 1712) — (0, 1), (0, 2) and (0, 3). 
These are the only states for which the minimum in the energy curves are less than the 
energy 2Ea = —323 eV of two isolated atoms in the ground state. Notice that as Z in- 
creases, the molecular electronic energy becomes larger than 2Ea, reflecting the fact that, in 
a superstrong magnetic field, forming such a tightly-bound molecule requires first activating 
one of the atoms to an excited state. However, as Z increases, the energy of the (mi, 7712) 
state does approach (mi )-!-£'„ (m2). Near the equilibrium separation, the energy obtained 
using HFMC agrees well with that of standard HF (the difference is less than 1%). Thus the 
standard HF is adequate for determining the equihbrium electronic energy of the molecule. 
However, the HFMC method is crucial to obtaining the correct large Z behavior of E{Z, 0), 
therefore the ahgned vibrational energy levels of the molecule (Sec. IV.B). 



These states can be similarly calculated using the HF method. Instead of Eqs. (3.9)- 
(3.10), the electron orbitals are $00 and $01 ■ Fig- 2 shows an example of the energy curve at 
B12 = 1- Clearly, the E[Z, 0) curve of such state is much shallower than those of the tightly- 
bound states discussed in Sec. HI. A. 2. In the limit of Z ^ 00, the energy curve approachs 
2Ea, i.e., no "activation energy" is needed to form a molecule in the weakly-bound state. 



E = ^{Hn + H22) - \[{H^i - H22f + 4.Hl^Y'\ 



(3.22) 



3. H2 Molecule: Weakly-Bound State {m,v) = (0,0), (0,1) 
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B. General Molecular Axis Orientation: R± ^ 



Unlike the case of Sec. III. A when the molecular axis coincides with the magnetic field 
direction, where we can obtain the interatomic potential E{Z, 0) with great accuracy, in 
the case when the molecular axis deviates from the magnetic field direction, the electronic 
energy E{Z^ R±) is much harder to calculate. This is because the azimuthal symmetry of 
the transverse wavefunction of an electron is broken. Although the electrons still stay in 
the ground Landau level, m in the Landau wavefunction W^m(r±) (Eq. [2.2]) is no longer a 
good quantum number, and the transverse wavefunction of an electron must involve mixing 
of different m-states. Nevertheless, we can still obtain a reasonable upper limit for the 
interatomic potential curve E{Z, R±), and hence an upper limit for the transverse vibrational 
excitation energy quanta huj±Q. We consider two ansatzs, appropriate for small R± and large 
R± respectively. 



1. Ansatz A 

Suppose the two protons are located at {±R±/2,0,±Z/2) in a rectangular coordinate 
system. For sufficiently small R±, the transverse wavefunction is expected to be close 
to Wmir^). Thus we assume the electron wavefunction in is given by $mo(r) = 
Wm{r±)fmo{z). The equation for fmo is the same as Eq. (3.4), except that the potential 
Vm{z) is replaced by 

/ Z R \ / Z R \ 

Vmm{z,R±/2) = Knmf k - — |, -^j + Vmm{\z + — |, -^j, (3.23) 

where 

V^m{z,R±/2)= fdh^lW^r^)]^- ^— - 

J |r — K_l/2| 

= J^^dq exp - q\z\y, (^^) Qg^^ , (3.24) 

(see Appendix A). Here Jq is the Bassel function of zeroth order and is the Laguerre 
polynormial of order m [^. We use a standard quadrature algorithm (e.g., [^) to evaluate 



Eq. (3.24). The Schrodinger equation similar to Eq. (3.4) can be solved to determine the 
eigenvalue emo{Z, R±), and the total electronic energy is then given by 

Emo{Z,R±) ^ emo(Z,Rj_) + j^^^-j-^^^^. (3.25) 

As noted before, in this general situation, m is not a good quantum number, but we never- 
theless use it to distinguish different electronic state. 

In this ansatz, the equations for H2 are also similar to those in Sec. III. A. We still assume 
the electron orbitals to be given by Eqs. (3.9)-(3.10). The HP equations (3.12)-(3.14) remain 
valid except the ion-electron interaction potential Vm{z) is replaced by Vjnm{z, R±/2). The 
electron-electron interaction kernels are unchanged. The total electronic energy is still given 
by Eq. (3.19) with e'^/Z replaced by eV(Z2 + i?^ )V2 

We now estimate the regime of vahdity of this ansatz. As an example, let us consider 
the ground electronic state of H2 . In general, the transverse wavefunction of the electron is 
a superposition of different Landau ground state wavef unctions, i.e., 

^±ir±)^^ArnWmM, (3.26) 

m 

and $(r) = ^±{r±)f{z) is the total wavefunction (see also Appendix B). For simphcity, just 
consider the first two terms in the expansion (3.26), i.e., $j_(r^) = AoWo{r±) + AiWi{r±), 
with I All \Ao\ for the ground state. Substitute $(r) into the Schrodinger equation and 
average over r^, we obtain (in atomic units) 

- - Voo{z)f + j-Voi{z)f = ef, (3.27) 

-^^/ - VnW + ^Vio{z)f = ef, (3.28) 

where V^^' is defined similar to Eq. (3.24). Since |Voi| <^ \Voq\ and |Voi| <^ l^iih from 
Eqs. (3.27)-(3.28) we have Ai/A^ ~ Vio/{Vu - Voo). Substitute this into the Eq. (3.27), we 

have 

-l^J-Voof + -^^^ f = ef. (3.29) 
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Comparing with the zeroth order eigenvalue (which does not take into account the 
mixing), the corrected eigenvalue for the ground state is then given by 

(0) , / VqiVw \ (0) , / VJjiKo \ /o o(\\ 

where (■ ■ ■) denotes expectation value. Requiring the second term to be smaller than the 
first, we have (Voi)^// <^ where we have used \e^Q \ ~ and \ef'' — ef^ ~ /. Since 

(Vol) ~-/Rx-V-\ r^-nj^) -R±p(\) --^R±, (3.31) 

\ r/oi \H/oi \H/ L^p 

the condition for the ansatz to be valid is R± <^ i^^'^p, i-e., the proton transverse displace- 
ment must be smaller than ~ p. 



2. Ansatz B 

At large the molecule should become two individual atoms (or atom plus ion). 
Here we set up a rectangular coordinate system so that the two protons are located at 
(0, 0, Z/2) and {R^_, 0, —Z/2). The electron wavefunction of 132" is assumed to be $^0(1") = 
W^m(r±)/mo(-2); i-G-, the clcctrou cloud is centered on one of the protons. Then the problem 
is essentially equivalent to calculating how an atom is affected by an external ion. The 
Schrodinger equation (3.4) still applies except that the potential Vm{z) is replaced by 

K^iz, R^) = V^(\z - + Vmm(\z + ||, (3.32) 

where the function Vm and Vmm are defined in Eqs. (3.6) and (3.24) respectively. The 
eigenvalue can again be solved and thus the total energy E{Z, R±) can be obtained. 

In this ansatz, since the electron wavefunction is not symmetric with respect to z —>■ —z, 
the numerical method used in Sec. III. A and III.B.l (see Paper I) need modification. Here 
we integrate the equation from both 00 and —00. The eigenvalue is obtained by matching 
the solution at 2; = (see |^). We also note that the classical quadrupole formula for the 
ion-atom interaction is not applicable here, since we always consider R± <, for a bound 
state. 
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For H2, we choose the two electron orbitals centered on each of the protons: 

^mA^)^Wm,iri.)fmAz), (3.33) 
$n^.(r) = Wm,{r^ - Rx)/m.o(;^)e-'^^-^/'. (3.34) 

The extra factor e^*^-'^-'-^/^ in ^^^(r) comes from a gauge transformation, so that the dis- 
placed Landau wavefunction (i"±— R-x) is still an eigenstate of the magnetic Hamiltonian 
with a fixed gauge (Eq. [3.2]), i.e., 

HB[WmAr± - Rx)e-^^-^/^x(i)] = 0. (3.35) 

With this ansatz for the basis wavef unctions, the HF equations given in Sec. III. A (Eq. [3.12]) 

can be applied, except that Vm{z) must be replaced by V^„^{z, R±) given in Eq. (3.32). Also, 
the direct and exchange kernels (Eqs. [3. 15]- [3. 16]) are replaced by 



(^1 - ^2, Ri.) = I d''r^^(fr2i_\W^,{vr^)\''\W^,{v2i. - Rx)!' — , (3.36) 

[zi - Z2, R±) = / rf'ri^rfWM^mi(rix)W^,(r2x - R±)W^;,(r2±)Vr;,(ri^ - R^) 

^^iBRAyi-y2)/2j__ /3_37N 

The function Drnim2 (^) R±) can be expressed as a sum of the function (see Appendix 
A) 

thus it can be evaluated using Eq. (3.24). For R± » p, the exchange interaction be- 
tween electrons can be neglected since the electron clouds are separated, i.e., we can set 
Eraim2{^i R±) — 0- Therefore, we only need to solve the corresponding Hartree equations: 



-K^m(2;, Rl) + —Km{z, R±) - Sr, 



where is given by 



fmoiz) = 0, m = mi, m2, (3.39) 



KmA^,R±) = ldz'U,o{zTDm,m2iz-z',R^), (3.40) 
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and similarly for ir^2- 

In Fig. 3, we show the energy curve for H2" at B12 = 1. The electron is assumed to be in 
the m = state. The electronic energy curves E{Z, Ro) are calculated using ansatz A with 
a fixed value of R± — Rq. Each curve has a minimum at Z — Z^q{R^. We see that this 
equilibrium position is almost independent of i?o, i.e., Zeq{Ro) — Zeq{0) = Zq. The curves 
E{Zo, R±) with a fixed value of Zg are calculated using both ansatzs discussed above. For 
R± less than a few times p, ansatz A yields lower energy, while for larger R±, ansatz B gives 
the correct behavior for the energy curve, i.e., E{Zo,R±) Ea as R_[_ increases. Similar 
behavior for H2 can also be obtained. It is evident from Fig. 2 that the curves E(^Zg, R±) are 
much steeper than E{Z, Ro) . Thus the molecule is tied much more "rigidly" to the magnetic 
field line than along the field axis. 

IV. RESULTS FOR THE MOLECULAR EXCITATION LEVELS 

In this section, we present our numerical results for the excitation levels of H2. The results 
for H2" are also included for completeness and for comparing with previous calculations (no 
previous results for H2 are available). 

A. Electronic Excitations 

The equilibrium electronic state is determined by the minimum in the energy curve 
E{Z,0) (cf. Fig. 1). For H^, the electronic state is characterized by a single quantum 
number m. For H2, there are two types of electronic excitations: the "tightly-bound" levels 
correspond to electrons in the (m, u) — (mi, 0) and (m2, 0) orbitals, and the "weakly-bound" 
excitation corresponds to (m, u) = (0, 0), (0, 1). We have calculated all the electronic bound 
states of and H2 for 0.1 < Bu < 10. The results for H2 are summarized in Table II 
(for the tightly-bound states) and Table III (for the weakly-bound state), while the results 
for B.2 are given in Table IV. Here, by "bound" we mean that the equilibrium electronic 
energy E^ of the molecule is less than E^ = -E'a(O), the energy of a single atom in the ground 
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state (for H^), or 2Ea, the energy of two atoms (for H2). Clearly, H2 has more electronic 
excitation levels than 132". As B increases, the number of bound levels in H2 increases. For 
B12 < 10, only single-excitation tightly-bound levels, i.e., those with mi = 0, are bound. 
The double-excitation levels, such as (1, 2) are not bound until the field strength increases 
to B12 ^ 50. Excluding the zero-point oscillation energy of the protons (see Sec. V), the 
dissociation energy the H2 molecule is given by D^°°^ = 2Ea — E^- 

We have also calculated the ground-state energy of the molecule in the stronger field 
regime. For Bu ^ 10, our numerical results can be well fitted to the form: 

Em ^ -0.091 (In 6)2-^ (a.u.). (4.1) 

More general fitting formula for D'^°°^ is given in Eq. (5.2). 

We note that as B increases, the energy \Em\ of the tightly-bound levels of H2 increases 
rapidly, while that of the weakly-bound level does not change appreciably. For Z = ln6 >> 1, 
the weakly-bound state is indeed an excited state of the H2 molecule. For B12 ^ 0.2, however, 
we find that the weakly-bound state actually has lower energy than the tightly-bound level 
(mi, 7712) = (0, 1). Thus for such relatively small magnetic field strength, the weakly-bound 
state is the actual ground state of the molecule. 

B. Aligned Vibrational Excitations 

In the standard Born-Oppenheimer approximation, the Hamiltonian describing the rela- 
tive motion of the protons is simply Hi = P^/ (2yu) + U {Z, R±), where /x is the reduced mass 
of the proton-pair, and the interatomic interaction potential U is given by the total electronic 
energy E{Z,R_i_), as calculated in Sec. III. For the vibrations along the 2;-direction, there 
is no magnetic force on the protons, and Eq. (4.1) is a good approximation. The aligned 
vibrations are governed by the potential U{Z,0) = E{Z,0), which we can fit to a Morse 
potential (e.g., p8| ) 

U{Z, 0) = Dm{l- exp[-P{Z - Z,)]Y + E^, (4.2) 
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where /3 is a fitting parameter, and 

D^ = U{oo,0)-E^. (4.3) 

Thus Djn — Ea{m) — Em for H2 , and D^. = Ea{mi) + Ea{m2) — Em for H2 (we consider the 
tightly-bound states only) . The ahgned vibrational energy levels are then given by 

E„,=/u.,(„,l + i)-^|^(„„ + i)', (4.4) 

where the vibrational energy quanta is 

nuj\\^h(3(^=^j . (4.5) 

The values of hui\\ and Dm for different bound electronic states and different magnetic 
field strength are given in Table II for H2 and in Table IV for Ha". In Fig. 1, the numerical 
E{Z,0) curve is compared with the fitted curve (Eq. [4.2]) for the (0,1) state of H2 at 
B12 = 1. We see that the fitting is indeed very good, especially for the bound region (below 
the dark line in Fig. 1). For H2, our results for hu!\^ are accurate to within about 5%. The 
Morse potential fits the E{Z, 0) curves of less well, but the resulting hu;\\ is still accurate 
to within about 10%. 

As discussed in Sec. II, ?iuj\\ (cf. Eq. [2.8]) is approximately proportional to (ln6)^//~^/^ 
times a slowly increasing function of B. Our numerical results confirm this approximate 
scaling relation. A better empirical scaling is ?iuj\\ oc (ln6)^/^/x~^/^. Thus for the (mi, 777.2) = 
(0, 1) state of H2 (/i ~ 918), we have 

hu;\\ ~ 0.13(ln6)^/V"^^^(a.u.) ~ 0.12(ln6)^/2 ^^y^^ ^fj^)_ ^4_g) 
For the ground state (m = 0) of we have 

hu\i ~ 0.085(ln6)5/2^-^/2^a.u.) ~ 0.076(ln6)^/2 ^^y^^ ^4 7) 

Both Eqs. (4.6) and (4.7) are accurate to within about 10%. These fitting expressions 
are indeed very satisfactory considering the approximation introduced when we use Morse 
potential to fit the numerical E{Z,0) curves. 
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There is no previous reliable calculation for H2 molecules. For H^, our results for the 
ground state electronic energy, interatomic spacing and aligned vibrational energy quanta 
fiojw agree with those obtained by Wunner, Herold and Ruder [|I^ , and those of Le Guillou 
and Zinn- Justin , who used a similar method as ours in the aligned cases. The slight dif- 
ference in huj\\ between our results and theirs is likely due to the different ways of extracting 
this quantity: we obtain it by fitting E{Z, 0) to a Morse potential, while they obtained it by 
evaluating the second derivative of E{Z^ 0) around the equilibrium separation. Le Guillou 
and Zinn- Justin also considered the effects of non-adiabaticity (i.e., mixing of different elec- 
tron Landau levels). This is negligible for field strength of interest in this paper (6 ^ 1). 



The variational calculation of Khersonskii |T6[ gave somewhat smaller (by about 20%) values 
for huj\\. This is due to the inaccuracy in his atomic binding energy. 



C. Transverse Vibrational Excitations 

Neglecting the magnetic forces on the protons, the transverse oscillations of the molecule 
are governed by the potential U{Zo,R±) = E{Zo,R±). Our calculation of this function is 
less accurate than the aligned case, and yields only an upper limit to the exact potential. 
For small-amplitude oscillation (see Sec. II), we fit this potential to a harmonic form 

6U{R^) = U{Zo,R±) - U{Zo,0) - i/i^ioi?l. (4.8) 

The transverse vibrational motion of the protons is therefore described by a two-dimensional 
harmonic oscillator. The numerical values for the transverse vibrational energy quanta are 
tabulated in Table II (for II2) and in Table IV (for H^). Only the results for the ground 
electronic states are given. 

Note that hujj_Q is larger than huJi^ for B12 ^ 1. Also, huj±o and huJi\ can be comparable 
or even larger than the electronic excitation energy spacings AE^- This is in contrast with 
the zero-field cases, where one has AEm >> humb >> Tiojrot- Although the actual values 
of huj±_o may be somewhat smaller than our results, the qualitative features revealed in our 
calculations are expected to be valid in general. 
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The discussion in Sec. II gives huj±Q ~ (^~^61n6)^/^/i~^/^ (Eq. [2.10]), where increases 
slowly with increasing B. Our numerical results confirm this approximate scaling relation 
and huj±o oc b^^'^{lnb)fi^^^'^ fits better the results in Table II and Table IV. For the (mi, 1112) = 
(0, 1) state of H2, we have 

/ h \'/' 

huj±o ~ 0.l25b^^^{\nb)fi-^/\a.n.) = 0.553 -— ln6(a.u.) = 0.11 b^/\\nb) (eV), (H2). 

\ bcrit J 

(4.9) 

For the m = state of H^, we have 

/ ft \'/' 

hiu±o ^ O.Ub^/^{\nb)fx~^^^{a.n.) = 0.62 -— ln6(a.u.) = 0.13 b^^^lnb) (eV), (H+). 

\ bcrit J 

(4.10) 

Our results for huj±o of also agree closely with those of Le Guillou and Zinn- Justin 



Tq obtained using their "static approximation", which is similar to ansatz A adopted in 
our paper (Sec. III.B.l). Their improved calculations indicate that the real value of huj±o 
can be lower by tens of percent (from about 10% for B12 = 0.1 to about 40% for B12 = 5). 
We expect our results for H2 to have similar accuracy. However, as noted in Sec. II. C, the 
present results apply only to the small-amplitude {R± ^ p) vibrations and relatively weak 
field {B <, Bcrit)- For sufficiently large field strength, the magnetic forces on the protons 
become important and can change the transverse vibration energy significantly, as we discuss 
below. 



V. EFFECTS OF FINITE PROTON MASS ON THE ELECTRONIC ENERGY 
AND MOLECULAR DISSOCIATION ENERGY 

Our calculations and results in the previous sections are based on Born-Oppenheimer 
approximation where the proton positions are fixed when we consider the electronic energy 
of the molecule. For finite proton mass, one can rescale the electronic energy by replacing the 
electron mass with an appropriate reduced mass. This only introduces a small correction (of 
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order me/mp), and is neglected in our paper. However, as noted in Sec. I, the separation of 
the proton and electron motion in strong magnetic field is much more complicated, especially 
in the regime of B ^ Bcru when the cyclotron energy of proton, 

hUp = h{eB/mpc) = (b/bcrit) In 6c«t (a.u.) = 6.3fii2 (eV), (5.1) 

is comparable with or larger than the spacing of the electronic energy levels. The "standard 
procedure" for separating the proton and electron motion leads to some ambiguities regard- 
ing the binding energy of H2 in the strong field regime; these are discussed in Sec. V.A. 
An alternative scheme, which is more suitable for B ^ Bcru, is described in Sec. V.B. An 
approximate expression for the "corrected" dissociation energy of H2 in the ground state is 
given by Eq. (5.7). 



A. Unbound States from the Standard Scheme When B >, B, 



crit 



In Sec. II and III we have followed the "standard procedure" for molecules, where one 
first considers the two protons as infinitely massive fixed at equilibrium separation Zq along 
the same field line, with their motion included only as an "afterthought". This is strictly 
valid only for S << Bcru, where the zero-point vibration amplitudes and the magnetic 
force on the protons are small. The two electrons, both in the lowest Landau level, are in 
cylindrical orbitals mi and m2 centered on the proton field line with radii given by equation 
(1.2). The Pauli principle requires mi 7^ m2 so that the ground state has mi + m2 = 1. 
As the p-p separation Z is allowed to increase, the system tends to two free H atoms, one 
in orbital state mi and another m2. The standard procedure for treating the two-body 



problem of hydrogen atom pO| , pl| deals with states where the transverse pseudomomentum 
of each atom is zero, in which case the protons must have Landau excitations mi and m2, 
respectively. The simplest state for the molecule with electronic orbitals mi and m2 is then 
the state where the protons have these Landau excitations at all separations (even though the 
transverse pseudomomentum is conserved only for the total molecule, not individual atoms). 
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This choice adds a Landau excitation energy (mi + m2)hu!p, i.e., even the electronic ground 
state would have an additional positive energy hup and would be unbound (relative to two 
ground-state atoms) for B » Bcru- This molecular state has the simplest wavefunction 
but not necessarily the lowest energy, since there are states where pseudomomentum is not 
zero (corresponding to finite separation of the guiding centers of the electron and proton; 
see Ref. [^). States without the additional energy hup are discussed in Sec. V.B. 

For infinite proton mass, the dissociation energy of H2 is given by D^°°^ = 2Ea{0) — Em- 
Our numerical results for the (mi, 7112) = (0, 1) ground state can be written in the following 
form 



where r varies slowly with 6 (r ~ 0.14 for b ~ 10'^ and r ~ 0.17 for b ~ 10^). For field 
strength b ~ bcrit the square bracket in Eq. (5.2) can be replaced by unity. As shown 
in Sec. II. B and Sec. IV. B, the aligned proton vibrations have an energy spacing of order 
huj\\ ~ /i^^/^Z}(°°) and a small vibration amplitude of order SZ ~ fi^^^^Zo, where Zg is the 
equilibrium separation between the protons (Eq. [2.5]). The inequality huj\\ « D^°°^ does 
not depend appreciably on the magnetic field strength, so for the ground molecular state 
we should be able to use hu!\\ in Eq. (4.6) for the aligned vibrations even when B » Bcru] 
furthermore, since 6Z << Zo, we do not need to consider the Pauli principle explicitly for 
the transverse wavefunctions of the protons. For treating this transverse motion, however, 
the magnetic force becomes important when b ^ bcrit, as can be seen from the ratio of 
the expressions in Eq. (5.1) and (4.9), hujp/huj±o — 1.81(6/6crit)"'^^^(ln6crit/ln&). We can 
give at least a plausibility argument for the inclusion into the "standard scheme" of the 
magnetic effects on the transverse motion from the following consideration: A free proton in 
the magnetic field B has a zero-point energy hujp/2. This can be thought of as the ground 
state energy in a "magnetic restoring potential" (l/2)mpu;p(/2_L/2)^, which gives a ground- 
state wavefunction of size R±_ ~ p (independent of mass) as in the Landau wavefunction 
(Eq. [2.2]). Thus the total restoring potential for for the protons in H2 is given by the sum 




0.2 



(5.2) 
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of the "electronic potential" SU{R±), which we have calculated in Sec. Ill, and twice (for 
two protons) of the magnetic restoring potential. For R± < p, we have seen that 6U{R±) 
can be approximated by the quadratic form in Eq. (4.8), so that the total potential is 
/^(^io + The size of the ground-state wavefunction is then <^ p, the approximation 

is justified and the excitation energy quanta is h{uj'j_Q + ujp)^^"^- Since the energy of hujp/2 
also exists in isolated H atom, the zero-point energy for the transverse oscillation of the 
molecule can be written as 

huj± = h{ujlo + ulY^'^ - hup. (5.3) 

For b « 

bcrit have hujj_ — hujj_Q as expected, and Eq. (5.2) shows that Tluj_\_I D^°^' ~ 
5.(6/6crit)^^^/ 1'^^ << 1- For h » bcrit, on the other hand, huj± ~ hujj_Q/{2ujp) ~ 0.016/^ 
so that huj±/D^°°'> ~ 0.15. Thus the transverse zero-point energy huj± remains less than 
the dissociation energy for the state given by the standard scheme, but the Landau energy 
hujp of the excited proton has to be added to the molecular energy also. The "corrected" 
dissociation energy in this scheme is then given by 

Clearly, Z}(**'^) becomes negative (i.e., the state is unbound) as b increases beyond bcrit- We 
shall see in Sec. V.B that an alternative scheme gives molecular bound states with lower 
energy for b > bcrit- 



B. The Alternative Scheme 



The alternative scheme we propose for the H2 molecule ground state is a generalization of 



the scheme for H atom described in Sec. IV of Ref. ||2^. In this scheme the transverse pseudo- 
momentum is not chosen as a good quantum number, and our approximate wavefunction 
will not be an exact eigenstate of the Hamiltonian. However, it does provide a suitable trial 
wavefunction and enable us to obtain a rigorous lower limit to the dissociation energy D. 
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In the absence of Coulomb interaction between the four particles, there are eight quantum 
numbers specifying the transverse degrees of freedom of the system: the Landau excitation 
number n and the orbital number m for each of the four particles. For the H2 ground state 
we then choose = for both electrons and both protons (so there is no hup contribution 
to the electronic energy of the molecule), and (mi, 7712) = (0,1) for the electrons. We 
can choose m = for both protons since, as mentioned, the proton 2;-wavefunction can 
be anti-symmetrized to satisfy the Pauli principle with little energy contribution. As a 
trial wavefunction we assume that the charge distribution of protons consists of two sheets 
separated in the 2;-axis by distance Z, with surface density given by |Wo(r±)p. Obviously, 
when the Coulomb potentials between the particles are restored, n and m for the individual 
particle cease to be good quantum numbers Q, but the trial wavefunction thus constructed 
will give an upper bound to the true ground-state energy of the molecule according to the 
variational principle. This "trial energy" can be calculated using the Hartree-Fock method 
described in Sec. II. A, subjected to two modifications: (i) The averaged electron-proton 
interaction potential Vm{z) in equation (3.5) is replaced by 

1 



Vm{z) ^ /f/Ve|W^m(r±e)r/c^%|W^o(r±p)p 



^ ( - f ) + + f 



where Dom(-2) is defined in Eq. (3.17); (ii) The proton-proton interaction term e^/Z in 
equation (3.19) is replaced by Dq^^Z) (although this modification has negligible effect on 
the energy except when Z ^ 0). 



^The only good quantum number for the transverse degrees of freedom is the total orbital angular 
momentum along the z-axis = sign(ej)(mj — n^), where sign(ej) = 1 for proton and sign(ej) = 



—1 for electron jS^]. For b » bcrit, the Landau excitation numbers n for both electrons and 
protons are "adiabatically" conserved and can be set to for the ground state. In this case 
Lz = nriip + m2p — rnie — m2e- Thus the true ground state of the molecule for B » Bcrit involves 
a mixing of many different (^T^le, ^T^2e, "m-ip, ^T^2p) states with the same L^. 
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The molecular energy -B^'*^ obtained by this alternative scheme is larger than the result 
Em obtained using the scheme of Sec. III. A (where the protons are treated as infinitely 
massive), by some amount Ae*^"'*^ The weakening of the electron-proton interaction is due 
to the spread of the proton wavefunction by an amount of order p. However, since Z » p 
the change involves only the logarithm of the Coulomb energy. This can be characterized by 
changing b in Eq. (4.1) to 6/(2C), where C is of order unity, i.e., ~ -0.091 [ln(6/2C)]2-^. 
To leading order in ln(2C), we then have 



as an alternative to Eq. (5.4). 

We have performed numerical calculations and found that the "trial" ground-state energy 
thus obtained agrees with the result using the scheme of Section 111. A to within 15%. For 
= 50, 100, 500, 10^ we found As^"'*) ~ 162, 191, 258, 294 eV, corresponding to C ~ 0.8 
for Bi2 — 50 and C ~ 0.9 for B12 — 10^. The numerical values for As^"'*) can be fitted by 
A£(°'*' ~ 0.06(ln6)2. This has the same scaling with b as hu± defined in Eq. (5.3). With this 
value of As'^"'*^ the dissociation energy given by Eq. (5.6) is larger than that from Eq. (5.4), 
and therefore represents the true molecular ground state for all b ^ bcrit- Our numerical 
calculation of Ae*^"'*) used a particularly simple trial wavefunction and a better wavefunction 
with the variational method would presumably lower Ae^"'*^ somewhat. This would lower 
the numerical value of b above which the true ground state is the state with no Landau 
excitations for either proton, obtained by the present "alternative scheme" . 

The fact that Ae^"'*) scales similarly with b as hu;± suggests that for practical purpose, 
the "corrected" dissociation energy of H2 in the ground state can be approximated by 



for all field regimes (& >> 1), where ?iuj± is given by Eq. (5.3). The numerical results for a 
wide range of field strength are summarized in Table V. 




(5.6) 




(5.7) 
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VI. CONCLUSIONS 



In this paper, we have, for the first time, studied and characterized the energy excitation 
levels of H2 molecule in a superstrong magnetic field {B ^ 10^^ G) which exists on the 
surfaces of many neutron stars. The main theoretical uncertainty of our calculations lies in 
the non-trivial separation of the motion of the protons and that of the electrons. Neverthe- 
less, we find that in such a strong magnetic field, H2 molecule exhibits completely different 
energy excitation levels as compared to its well-known zero-field counterpart. The fact that 
the excitation energies associated with the oscillations of the protons are comparable to the 
electronic excitations indicates that the statistical weight of a H2 molecule is not much larger 
than that of a H atom. This greatly simplifies the calculations of the chemical equilibria of 



various forms of H in a neutron star atmosphere [IS 



Larger hydrogen molecules and chains can also form in a superstrong magnetic field. 
Their ground state binding energies have been calculated in Paper I. It is expected that 
these larger molecules possess qualitatively similar energy excitation levels as those of H2 
considered in this paper, with one exception: For a long chain molecule H„ with 1 << 
n « [6/(ln6)^]^/^, the spacing Zo along a field line between adjacent protons decreases 
with increasing n approximately as n~'^. The fractional zero-point vibration amplitude 
AZ/Zo is of order (m^./mpY^^n^^'^. The aligned vibrations thus become more pronounced 
as n increases (and can lead to "internal pycnonuclear reactions" which will be discussed in 

H). 

There is no question that the exotic molecules considered in this paper exist on the 



surfaces of some neutron stars with B12 ^ 10^^ G and temperature T ~ 10^ — 10^ K [jl8 
For very low surface temperature (T < 10^ K), the atmosphere is likely to condensate 
into a metallic state, since the hydrogen metal has the largest binding energy. However, 
for the astrophysically more interesting temperature range (T ^ 10^ K), the outer layer 
of a neutron star will predominantly exists in the form of nondegenerate gas of individual 
atoms and small molecules: e.g., when T ~ 3 x 10^ K, the photosphere of a neutron star 
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is dominated by atoms if B12 = 1, while it is dominated by H2 if B12 = 10. The existence 
of H2 in the atmosphere will give rise to appreciable radiative opacity. For example, since 
the proton separation in H2 is different from that in (see Table II and Table IV), the 
photo-ionization cross-section from the ground state of H2 is expected to be small according 
to the Franck-Condon principle. However, photo-ionization from an excited vibrational state 
or electronic state, for which the proton separation is close to that in Hg" ground state, can 
provide significant continuum opacity. These issues may warrant further study, especially 
in light of the increasing possibility of the spectroscopic studies of isolated neutron stars by 
future X-ray/EUV satellites. 
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Fellowship in theoretical astrophysics at Caltech. 

APPENDIX A: COULOMB INTEGRALS FOR LANDAU WAVEFUNCTION 

In this Appendix, we derive Eqs. (3.24) and (3.38). First consider the function 
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m). 



(Al) 



r — r, 



o 



with r, 



o 



ToX. Since 




(A2) 



we have 




(A3) 



Using the general result for the matrix element [| 
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{m! > m), (A4) 

where 9g specifies the angle of in the — Qy plane, and is the Laguerre polynomial 
of order n we have 



) = e-^y'Lj^l). (A5) 



Substitute Eq. (A5) into (A3), and integrate out dqz and d9q using 



and 



J dOge-''^^'^''^"'^'' = 27iJo{q^ro), (AT) 



we obtain 



i.e., Eq. (3.24). Note that using Eq. (A4), a more general expression can be obtained for the 
matrix element 

Vm'm{z,ro) = {m'\- — ^ — -\m) 



Next consider Drnim2 defined in Eq. (3.36). Changing variable (r2_L — R_l) r2±, we 
have 

Dm,mM-Z2,R±) = fdhi^dh2±\WmAri^)\^\Wm,{r2±)\^-. ^ ^ , . (AlO) 

J |ri — r2 — -rt_L | 

Using Eq. (A2), we have 

Dm,m,{z,R±) = 7^ /^e-*'i-^^e*^^-(mi|e^'i^-'-^|mi)(m2|e*'i^-'-^|m2). (All) 



34 



Again, using Eq. (A4), and integrating out dqz and dOq with Eqs. (A6)-(A7), we obtain 

Dm,m,{z,R±) = /^"dge-^'-<'I^IJo(gi?±)L„,(^g')L^,(^g'). (A12) 
Now defining the coefficient (mi, 777,2) via (see Paper I) 

Lrm{-^)Lm2{-^) = ds{mi,m2)Ls{x), (A13) 

Eq. (A12) then becomes 

mi+m2 



which reduces to Eq. (3.38) after using Eq. (A8). 



APPENDIX B: MORE ACCURATE CALCULATION OF H 



2 



An "exact" treatment of H2' for general orientation of the molecular axis proceed as 
follows. Consider the coordinate system of Ansatz A in Sec. III.B.l. When b » 1, the 

most general electron wavefunction for the 1/ = state can be written as 

^mo{r)=J2^m{r±)fmo{z). (Bl) 

m 

Substituting this into the Schrodinger equation and averaging over the transverse direction, 
we obtain a set of differential equations for fmo{z): 

fmo{z)-—Y,^rnm'{z,R±/2)fm'o{z)=S^ofmo{z), m = 0, 1, • • • (B2) 



2meP^ dz^ p 
where Vmm' is defined similar to Eq. (3.23): 

and the function Vmm' can be evaluated using Eq. (A9). Equation (B2) is subject to the 

boundary conditions dfmo/dz = at 2; = and fmo — > as 2; — >• 00. The normalization 
condition requires 

35 



E/ dz\Uiz)\'-i- (B4) 

The set of equations (B2) can be solved numerically using an iterative scheme similar to 
that used for solving the Hartree-Fock equation (Paper I). Successively accurate results can 
be obtained by using increasing number of terms in the sum in Eq. (Bl). The lowest energy 
state corresponds to the solution satisfying 

/oo roo roo 

dz\fooiz)\'> dz\hoiz)\'> dz\f2oiz)\'>---. (B5) 
-oo J—oo J —oo 

Generalization of this method to H2 molecule is much more complicated. 
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FIGURES 

FIG. 1. The electronic energy curves E{Z,0) for the tightly-bound states of H2 molecule at 
B = 10^^ G when the molecular axis is aligned with the magnetic field axis. The electrons 
occupy the {m,v) = (mi,0) and (m2,0) orbitals (mi 7^ 777,2). The solid line is for the state 
(mi, 7712) = (0,1), the short-dashed line for (0,2), the long-dashed line for (0,3). The dotted line 
is from the fitting using the Morse potential (Eq. [4.2]). The dark solid line corresponds to the 
energy of two isolated H atoms in the ground state 2Ea = —323 eV. 

FIG. 2. The electronic energy curves E{Z, 0) of H2 molecule at S = 10^^ G when the molecular 
axis is aligned with the magnetic field axis. The solid line corresponds to the tightly-bound state 
in which the electrons occupy the {m,u) = (0,0) and (1,0) orbitals, the dashed line corresponds 
to the weakly-bound state in which the electrons occupy the (0, 0) and (0, 1) orbitals. 

FIG. 3. The electronic energy curves for the ground state of at S = 10^^ G. The light 
lines show the E{Z, Rg) curves with a fixed R± = Rp for i?o = (solid line), Ro = p (dotted line) 
and Ro = 2/5 (dashed line). The dark lines show the function E(Zo,R±) for a fixed value of Zq 
given by the equilibrium separation of the protons. The solid line is calculated using ansatz A, the 
dotted line using ansatz B (see Sec. III.B). 
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TABLES 

TABLE I. Energy levels Ea{m) (in eV) of hydrogen atom in superstrong magnetic field. 
The levels are specified by the quantum number m, while the longitudinal node v = {). Here 

5i2 = 5/(10^2 Q) 



B\2 










^a(4) 




0.1 


-76.4 


-52.5 


-43.3 


-38.0 


-34.4 


-31.8 


0.5 


-130.2 


-92.8 


-77.8 


-69.0 


-63.0 


-58.6 


1 


-161.5 


-116.9 


-98.7 


-88.0 


-80.6 


-75.1 


2 


-198.5 


-145.8 


-124.1 


-111.2 


-102.2 


-95.5 


5 


-257.1 


-192.6 


-165.5 


-149.2 


-137.8 


-129.2 


10 


-309.6 


-235.1 


-203.5 


-184.3 


-170.9 


-160.7 
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TABLE II. The tightly-bound energy levels of H2 molecule in which the electrons occupy the 
(m, I/) = (mi,0) and (m2,0) orbitals (mi 7^ 7712). Here B12 = B/{10^^ G), Ea is the ground state 
energy of H atom, mi,m2 are the quantum numbers specifying the electronic excitations, is 
the electronic energy of the molecule, Zg is the equilibrium interatomic separation (oq is the Bohr 
radius), Dm is defined by Dm = Ea{mi) + Ea{m2) — Emimi,m2), huj^^ is the aligned vibrational 



energy quanta, and huj±o is 


the transverse vibrational 


energy quanta 


(neglecting 


the magnetic 


forces on protons). 














Bi2 2Ea (eV) 


mi, 7712 


Em (eV) 




Dm (eV) 


nu\\ (eV) 


hu±o (eV) 


0.1 -152.8 


0,1 


-161 


0.52 


31.7 


3.0 


2.6 


0.5 -260.4 


0,1 


-291 


0.30 


67.5 


7.2 


8.7 




0,2 


-264 


0.32 


55.7 


6.3 




1 -323.0 


0,1 


-369 


0.25 


91.0 


9.8 


14 




0,2 


-337 


0.26 


76.5 


8.8 






0,3 


-323 


0.26 


73.0 


8.3 




2 -397.0 


0,1 


-466 


0.20 


121 


13 


23 




0,2 


-425 


0.21 


103 


12 






0,3 


-408 


0.21 


98.3 


11 






0,4 


-398 


0.22 


96.8 


11 




5 -514.2 


0,1 


-623 


0.15 


173 


19 


42 




0,2 


-573 


0.16 


150 


18 






0,3 


-550 


0.16 


143 


17 






0,4 


-537 


0.16 


142 


16 






0,5 


-527 


0.17 


141 


16 






0,6 


-519 


0.17 


140 


16 




10 -619.2 


0,1 


-769 


0.12 


224 


25 


65 




0,2 


-709 


0.13 


196 


23 






0,3 


-682 


0.13 


188 


22 





42 



0,4 


-666 


0.14 


185 


22 


0,5 


-654 


0.14 


183 


21 


0,6 


-645 


0.14 


183 


21 


0,7 


-638 


0.14 


182 


21 


0,8 


-632 


0.14 


182 


20 


0,9 


-627 


0.14 


182 


20 


0,10 


-623 


0.14 


181 


20 


TABLE III. The energy of the weakly-bound state of H2 in which the electrons occupy the 


(m, u) = (0, 0) and (0, 1) orbitals. 


Here Bi- 




G), Em is the energy of the molecule. 


= \Em\ — 2\Ea\ is the dissociation enerj 


5y of the level (neglecting the zero-point oscillation 


energy of the protons), and Zq is the equihbrium interatomic separation (ao 


is the Bohr radius). 


B12 E^ (eV) 






(eV) 


Zo (flo) 


0.1 -167 






14 


1.5 


0.5 -279 






19 


1.3 


1 -344 






21 


1.1 


2 -421 






24 


0.99 


5 -542 






28 


0.89 


10 -649 






30 


0.72 
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TABLE IV. Bound-state energy levels of H+. Here B12 = 5/(10^^ q)^ jg ^he ground state 
energy of H atom, m is the quantum number specifying the electronic excitations of the molecule, 
Em is the electronic energy of the molecule, is the equilibrium interatomic separation (oq is the 
Bohr radius), Dm is defined by Dm = Ea{m) — Em, h^^w is the aligned vibrational energy quanta, 
and fiu±_Q is the transverse vibrational energy quanta (neglecting the magnetic forces on protons). 



-B12 


Ea (eV) 


m 


Em (eV) 




Dm (eV) 


huj\\ (eV) 


^u;_Lo (eV) 


0.1 


-76.4 





-99.9 


0.62 


23.5 


2.0 


3.1 


0.5 


-130.2 





-182 


0.35 


51.8 


4.9 


9.8 


1 


-161.5 





-232 


0.28 


70.5 


6.6 


16 






1 


-162 


0.40 


44.8 


4.4 




2 


-198.5 





-293 


0.23 


94.6 


9.0 


25 






1 


-207 


0.32 


61.5 


5.9 




5 


-257.1 





-393 


0.18 


136 


13 


45 






1 


-284 


0.24 


91.1 


8.6 




10 


-309.6 





-486 


0.15 


176 


17 


70 






1 


-356 


0.19 


121 


12 
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TABLE V. The dissociation energy of H2 molecule in the ground state in a superstrong mag- 
netic field. is the dissociation energy assuming infinite proton mass, while D includes the 
(approximate) correction of the molecular zero-point energy (Eq. [5.7]). B12 = i3/(10^^ G), h(jj\\/2 
is the zero-point energy for the aligned vibration, and hu± is the zero-point energy for the transverse 
vibration (Eq. [5.3]). Note that for B12 < 0.2, the ground state is actually the "weakly-bound" 
state (see Sec. 11. A and IV. A), and the zero-point energy has been neglected. All energies are 
expressed in eV. 



B12 


0.1 


0.5 


1 


5 


10 


100 


500 




14 


31 


46 


109 


150 


378 


615 


D 


~ 14 


21 


32 


79 


110 


311 


523 


^a;||/2 




3.6 


4.9 


9.5 


12 


22 


31 


Tiu± 




6.1 


9.1 


21 


28 


45 


61 



45 




0.2 0.4 0.6 0.8 



